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Abstract 

In this paper, the Hamiltonian formulations along with the Poisson brackets for two- 
dimensional (2D) electron magnetohydrodynamics (EMHD) flows are developed. These for¬ 
mulations are used to deduce the Beltrami states for 2D EMHD flows. In the massless 
electron limit, the EMHD Beltrami states reduce to the force-free states, though there is no 
force-free Beltrami state in the general EMHD case. 
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1. Introduction 


In the magnetohydrodynamics (MHD) model, the dynamics is dominated by ions with 
electrons serving to shield out rapidly any charge imbalances. In electron MHD (EMHD), 
with i ^ ps^, s = i,e, ps being the gyro-radius, on the other hand, the dynamics is dominated 
by electrons with the demagnetized ions serving to provide the neutralizing static background 
(Kingsep et ah |T], Gordeev et ah 0 ). The assumptions underlying the EMHD model are 
i -C di, where dg = c/up^ is the skin depth, and that the frequencies involved are greater 
than cjci and O 0 p^, Uc being the cyclotron frequency. Some theoretical developments including 
two invariants for the 2D EMHD equations were given by Biskamp et ah [5] . 

On the other hand, equations governing several plasma dynamics models are known to 
admit a significant class of exact solutions (Shivamoggi [1]) under the Beltrami condition - 
the local current density is proportional to the magnetic field - the force-free state (Lundquist 
[5], Lust and Schluter [6]). These Beltrami solutions are known to correlate well with real 
plasma behavior (Priest and Forbes [7], Schindler [8]). The purpose of this paper is to 
develop the Hamiltonian formulation along with the Poisson bracket for two-dimensional 
(2D) EMHD flows. This formulation is used to deduce the Beltrami states for 2D EMHD 
flows. These considerations are extended to the massless-electron limit (Huba and Rudakov 
[S], Rheinhardt and Geppert [10], Wood et al. in])- 

2. Governing Equations of 2D EMHD 

On assuming that the displacement current dE/dt is negligible, which is valid if 
u! <C Up^/oJce, the continuty of the electron flow implies 

^ = -V ■ (neVe) = V ■ (-) = - V ■ (V X B) = 0 
at \e J e 

or 


rie = const. (1) 

The 2D EMHD system of equations can then be written in terms of two scalar potentials 
(Biskamp et al. [3]) - the magnetic flux function A describing the in-plane magnetic field 
B = V X A iz and the stream function describing the in-plane solenoidal electron flow 
velocity Ve = V x -0 i^, which is proportional to the in-plane current density (so also 
represents the out-of-plane magnetic field): 

* the equation of generalized vorticity: 

|(, + ^^)+(,,.V).-^(B.V)J = 0 (2) 

* the generalized Ohm’s law: 

k(yl + ^j) +(v..V)(2l + ij) =0 (3) 


where. 
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( 4 ) 


-J = -VM , a; = -VV- 
c 

Equations (1) and (2) may be rewritten as 


where, 


d'lpe 

dt 




BA, 

dt 


[■0, ^e] 


( 5 ) 

( 6 ) 


V'e = rneTle {i>/ dg + Cu) '| 

A, = A+{dl/c^)j [ (7) 

[A, B] = VAxVB ■ 1,. J 

In the linear limit, on assuming a constant magnetic held Bq in the x-direction and a 
space-time dependence for the perturbed quantities, equations (5) and (6) give the 

whistler waves [i/jk = ^kA^ - the whistler condition). 


u ^ k‘^deVA, ( 8 ) 

which are dispersive unlike the Alfven waves. Here, Va^ = Bq/ y/m^n^. 


3. The Hamiltonian Formulation 


Equations (5) and (6) have the Hamiltonian formulation (Biskamp et ah |3]), 


H = 


1 

2 


'0'0e H— Af,J ) dA 


( 9 ) 


A 

A being the area occupied by the plasma. 

The Hamilton’s equations are, on taking {'il>e,Ae) to be the canonical variables, given by 


/ d^/je\ 


t6H\ 

dt 

= 3 

6lpe 

dA, 


6H 

\ dt J 


\5aJ 


( 10 ) 


where ^ is a {‘ipe, He)-dependent differential operator which produces a skew-symmetric trans¬ 
formation of vector functions vanishing on dA and satisfies a closure condition on an asso¬ 
ciated symplectic two-form. 


/[■■V'J KAlX 

VhA] 0 / 

and 6H/Sq is the variational derivative. Equation (10) is just equations (5) and (6). 


( 11 ) 
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The operator 3 induces a Poisson bracket as follows, 


[C G] ^ 


A 


A 


6F 6F 

SAe 

6F 6F 
5lpe ’ 

5F 5G 


, a 

/ 


V 


5G 5G 

6G 


6ijj, 


-,'06 


+ 


5G 


' 5G , 

Jaa^^ 


,^6 


0 


_(50e’ 506. 


+ ^6 


506 
' 5F ^ 
50e’ 5Ae 


5G 5F 


( 12 ) 


_50e’ 5Ae_ 




which is a bilinear function defined on admissible functionals F [0e, Ae\ and G [0e, A^ satis¬ 
fying the gauge condition, 


V- 

and the boundary conditions, 

6F 


506 


/6F\ 


/6G\ 

506 

= 0, V- 

506 

6F 


SG 

\6Ae/ 


\6aJ 

6F 

1 r 5G 

5G 


= 0 in ^ 


(13a) 


5216 


506 


5216 


= 0 on dA. 


(13b) 

(14a) 

[[F, G],K] + [|G. K], F] + [[K, F]. G] = 0. (14b) 

In the absence of elecron inertia (dg 0), (12) reduces to the MHD result given by Holm 
et ah 


The Poisson bracket (12) satisfies the anti-symmetry property, 

[F,G1 = -[G,F] 

and the Jacobi identity. 


4. Beltrami States 


The Casimir invariants for this problem are annihilators (with respect to any pairing 
functionals) of the Poisson brackets, which become degenerate when expressed in terms of 
these natural quantities. The Casimir invariants are therefore solutions of the equations. 


a 


506 

56 

\IaJ 



(15) 
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It may be verified that two such solutions are, 


6 ( 1 ) = 

[ f (Tie) dA 

(16) 

A 


6 ( 2 ) = / 

rjj^g (He) dA 

(17) 


A 

for any functions / and g. (16) and (17) were deduced by Biskamp et ah [3] by working 
directly with equations (5) and ( 6 ). C(i) signifies the conservation of generalized magnetic 
flux in EMHD. 

Minimization of if, keeping C(i) fixed, gives 


SH _ 56(1) 

Me Me 


(18) 


1J = A,i,/'(Ae) (19) 

which prescribes the Grad-Shafranov equilibria for 2D EMHD. There is no force-free state 
in 2D EMHD. This situation changes, however, in the massless-electron limit (see Section 
5). 

On the other hand, minimization of H, keeping 6 ( 2 ) fixed, gives 


SH _ 56(2) 

6i)e Si’s 

or 


( 20 ) 


'P = ^I2)SW 

or 


( 21 ) 


where. 


Ve = X(2)g' {Ae) Be = oBe 


Be = V X He iz = B — dgV^B. 
(21) describes the generalized Alfvenic state in 2D EMHD. 


( 22 ) 

(23) 


5. EMHD in the Massless Electron Limit 

In the massless electron limit, the Ohm’s law is 

Eh - VeXB = 0 

c 


(24) 
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which corresponds to the constraint of zero electric field in the reference frame moving with 
the electron flnid. 

Using, 


Ve =-= —X B (25) 

rieC UeC 

and Faraday’s Law, 

-if = V><E (26) 

equation (24) leads to (Huba and Rudakov [9]), 

r 

+-V X [(V X B) X B] = 0 (27) 

dt UeC 
or 


/\ • c* 

+ —B, {d,B,) = —B^ (a,R,). 
dt n.e n^e 


(28) 


In 2E0, where A = Ai^, equation (28) becomes (upon normalizing t appropriately). 


-^ = [B..A]. (29) 

Equation (29) implies that the out-of-plane magnetic field B^ plays the role of the stream 
function for the in-plane electron flow, as to be expected. 

Equation (29) has the Hamiltonian formulation. 


Ff 


JABz dA. 

A 


The Hamilton’s equation, on taking A to be the canonical variable, is given by 


(30) 


dA 

Ja 


(31) 


where 3 is an A-dependent differential operator which produces a skew-symmetric transfor¬ 
mation of scalar functions vanishing on dA and satisfies a closure condition on an associated 
symplectic two-form. 


3=l;A]^ 


Equation (31) is then just equation (29). 

The operator 3 induces a Poisson bracket as follows. 


(32) 


^See Appendix for the 3D formulation in the massless-electron limit. 
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A 


6F 

Ja 


5G 

Ja 




dA 


( 33 ) 


A 


A 


Ja’m 


dA 


which is a bilinear function defined on admissible functionals F[A\ and G[A], satisfying the 
conditions, 


V ■ = 0 in A and |^|, |^| = 0 on dA. 

\dG/6AJ ' 6A' ' 5A' 

The Casimir invariants are then solutions of the equation, 


(34) 



(35) 

= j f(A)dA 

(36) 


which leads to 


f{A) being an arbitrary function. 

Minimization of TC, keeping C fixed, leads to the Beltrami state, 

or 

B,= 

which leads to. 


describing a force-free state in the massless-electron limit (contrary to the general EMHD 
case). So, EMHD, in the massless-electron limit, mimics usual MHD (this result holds in 3D 
as well - see Appendix). 



(37) 

Af (A) 

(38) 

(A) B = aB 

(39) 


6. Discussion 

The intrinsic tendency of real plasmas (or fluids) toward Beltramization is a long standing 
issue awaiting full clarification. Work to date has indicated that the Beltramization process 
provides the means via which plasma systems can accomplish. 
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* an effective depletion of non-linearities, 

* ergodicity of the streamlines of the plasma flow in qnestion (Moffatt [T3]). 

* selective dissipation of total energy (Woltjer [S]). 

In this paper, the Hamiltonian formnlations along with the Poisson brackets for 2D EMHD 
flows are developed. These formulations are used to deduce the Beltrami states for 2D EMHD 
flows. Though there is no force-free Beltrami state in 2D EMHD flows, in the massless- 
electron limit Beltrami states become force-free. So, EMHD, in the massless-electron limit, 
mimics usual MHD (this result holds in 3D as well - see Appendix). 
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Appendix 

In the massless electron limit, upon appropriately renormalizing t, equation (27) can be 
rewritten as 


aB 


+ V X {(V X B) X B} 


Equation (AT) has the Hamiltonian characterization. 


0 


(AT) 


= 



(A-2) 


where V is the volume occupied by the plasma. Hamilton’s equation, on taking B to be the 
canonical variable, is given by 


dB _ Adf 


(A-3) 


where ^ is a B-dependent differential operator, which produces a skew-symmetric transfor¬ 
mation of vector functions vanishing on dV and satishes a closure condition on an associated 
symplectic two- formll. 


a = Vx {Bx (Vx (■))}. 


(A-4) 


^(A-4) is similar to the one given by Olver m for incompressible hydrodynamics (with the flow vorticity co taking the place 
of B). 



Equation (A-3) is then just equation (A-1). 


The operator S may be seen to induce the Poisson bracket (similar to the one given by 
Kuznetsov and Mikhailov jl6j). 




V X 


6G 


X V X 


6F 


dV 


(A-5) 


which is a bilinear operation defined on the admissible functionals satisfying the conditions, 


f6F/6B\ 
■ [^5G/6B J 


= 0 in V 


and 


SF 6G 


0 on dV. 


The Casimir invariant is then given by the solution of the equation, 


which leads to the magnetic helicity. 



= 0 


e 


j A-BdV. 


V 


Minimization of TC, keeping C fixed, leads to the Beltrami state. 


m ^ 
^ ~ 


(A-6) 


(A-7) 


(A-8) 


(A-9) 


B = AA (A-10) 

which leads to 

J = cAB (ATI) 

describing again a force-free state in the massless-electron limit. So, EMHD, in the massless- 
electron limit, mimics usual MHD. 
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